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We investigate the influence of electron-electron interactions on the density of states of a ballistic 
two-dimensional electron gas. The density of states is determined nonperturbatively by means of 
path integral techniques allowing for reliable results near the Fermi surface, where perturbation 
theory breaks down. We find that the density of states is suppressed at the Fermi level to a finite 
value. This suppression factor grows with decreasing electron density and is weakened by the 
presence of gates. 
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The suppression of electron tunneling into a conductor 
at low bias voltages, a phenomenon known as zero-bias 
anomaly (ZBA), has been under consideration theoret- 
ically and experimentally for more than 20 years [1-8]. 
This effect can be related to a reduction of the (tunnel- 
ing) density of states (DOS) at the Fermi level, which is 
a clear signature of interaction effects that are otherwise 
often disguised in good conductors. 

In the past the focus has been mainly on disordered 
systems with slow diffusive electron motion enhancing 
Coulomb effects so that strong ZBAs arise 01- More re- 
cently^attention has shifted to one-dimensional ballistic 
wires '9| , where the Tomonaga-Luttinger model leads to 
a power law suppression of the DOS [lOj. Similar to dis- 
order, the reduction of dimension hinders fast spreading 
and forces the particles to interact more strongly. Two- 
dimensional ballistic systems are a borderline case be- 
tween strongly interacting one-dimensional fermions and 
weakly interacting three-dimensional systems. The DOS 
is also of special interest, since two-dimensional electron 
gases (2DEGs) arise in a variety of semiconductor mi- 
crostructures, such as GaAs heterostructures and Si in- 
version layers that have attracted a lot of attention lately. 
Further, 2DEGs provide electrodes in experimental tun- 
neling setups measuring properties of ID quantum wires. 
Since such experiments do not test the DOS of the wire 
alone, it is important to know the form of the DOS of 
the electrode. 

In earlier work [HI , based on perturbation theory 
in the Coulomb interaction, the change of the DOS of 
2DEGs in the absence of a screening gate was found to 
show a cusp Sv{€)/vo = |e|/4eF at the Fermi edge. This 
prediction is rather irritating since it is independent of 
the strength of the Coulomb interaction. Here we reex- 
amine the problem employing nonperturbative path inte- 
gral techniques and show that the suppression of the DOS 
at the Fermi level sensitively depends on the electron den- 
sity and screening by gates. In contrast to the ID case, 
in the absence of disorder, the DOS at the Fermi level 
remains finite also at zero temperature yet reduced from 
the bare density that is approached at larger energies. 



The suppression increases with the effective interaction 
strength, i.e., with decreasing electron density. 

The interacting electron gas can be described by the 
action for a fermion field coupled to an electric poten- 
tial (j) 

+ e/|f,|^^(p)0(p-p')^(p')} 

+ / ^V;(p)G„-i(p)^(p), 

where Vq is the Coulomb interaction potential and the 
operator GQ^{p,t) — idt ~ e(p) (we set h ~ I) contains 
the electronic dispersion relation e(p) — /2m* with the 
effective mass m*. The time integration path is along 
the Keldysh contour. For simplicity we restrict ourselves 
to spinless Fermions but account for spin degeneracy by 
appropriate factors of 2. Further, we add a source term 
Jtp + i/jJ to the action allowing us to calculate the Green 
function later. 

The Fermion fields can be integrated out in the stan- 
dard way |l3| yielding an effective electromagnetic action 

5eff[0] = l^tf^-Y J </'(-q)V'(q)</'(q) 

+Trln [G,\p)S{p - p') + e0(p - p')] 
+ / J(p)GM(p)J(p)| , (1) 

where the trace means Tr — J ^^J^2 (27^)2 S{p — p') and 
G[0] = [Goi+e0]-i. 

Expanding the logarithm in the action with respect to 
the Coulomb field to quadratic order, the first two terms 
in Eq. |^ can be combined to the Gaussian action of the 
electric potential field 

Sf^ = -y^^*/ (02 H-<Dv-\cimci) , (2) 
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where 



(3) 



is the dynamicaUy screened interaction. Separating fields 
on the upper and lower parts of the Keldysh contour by 
introducing doublets, we can pass over to a convenient 
representation in terms of retarded and advanced func- 
tions (and mixtures thereof). Then Pq becomes a matrix 
in Keldysh space, with a retarded part given by 



d'^P "-F(£(p + q)) -nF(£(p)) 
(27r)2 uj + i7]- e(p + q) + e(p) ' 



(4) 



where rtF(e) = [1 + e^*^]^^ is the Fermi distribution func- 
tion at inverse temperature /?. The Gaussian approx- 
imation made in deriving Eq. ^ is equivalent to the 
random phase approximation (RPA), which is well es- 
tablished for high electron densities Hs corresponding to 
small values of the Brueckner parameter < 1, but usu- 
ally gives still very reasonable results for larger r^. As 
standard, Ug = l/7r(rsao)^ with the effective Bohr radius 
flo — £d/'7i*e^, where Sd is the dielectric constant. 
The quantity of interest here is the DOS 
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ImG^(p,e) 



(5) 



where is the retarded Green function. The Keldysh 
matrix Green function is obtained as the mixed second 
order functional derivative of the partition function with 
respect to the sources yielding 



G{r~r',t-t') ^ / V<l>G[cl)]{r-r',t~t') e 



Following Schwinger [isl we try to find a func- 
tional k[(j3]{x,t) describing a local gauge transformation 
'ipix,t) e*'=['^I(^'*)V(a;,i) such that 

G[(l)]{x-x\t-t') 

= e*'=['*l(^'*)Go(.T ~x',t- t')e-*M0](^',t') . (g) 

Linearizing the dispersion relation near the Fermi sur- 
face, we can write for w <C ep 

[idt - ep + iv • V + e(/)(a;, t)]G[(/)] {x - x\t ~ t') 
^S{x~x\t~t') (7) 

where v = p/m* is the Fermi velocity. Eqs. © and O 
determine A: as a linear functional of 4>. For the Green 
function G-^ — —i{ip{x,t)'ip{x',t')) we then obtain 

G>{x-x',t-t')^G^{x~x',t- t')e-^(^-^' 

where Gq is the free Green function and J{x — x' ,t — t') 
is determined by the remaining Gaussian path integral 



over the (p fields. This function, which we need for equal 
space arguments only, may be written as 



J{t) = J{x = x',t) 
where 

y(^) = - 



duj 
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ImF(w) , (8) 



V\q,L0). (9) 



(27r)2 {lo + 11] — V ■ q)2 
Now, from Eq. we finally obtain the formal result 



^(e) =1^0 de' 



-fie 



1 + e-l^'^' 



PU 



(10) 



where — m/i: is the bare DOS and we introduced the 
spectral density 



P{e) 



dt 
2^ 



(11) 



To evaluate the DOS explicitly, we first note that the 
polarization function Q reads for small q 



27r 



d6 vq cos 6 



(12) 



P^{(l,uj) = -1^0 I ^ 

Jq ztt S2 — vqcosb 

where fl — uj + irj. If we introduce the function 
giq) = [l + iyoVoiq)]-'/\ 

which is unity in the noninteracting case and less than 
one in presence of interactions, Eq. ((SJ may be combined 
with Eqs. Q and ((T^ to read 



2,71 



dq 



nqg\q)Vo{q) 



W-{vqn [n + vq)J^~[l~g^{q)]n 



The first factor in the denominator represents the parti- 
cle pole and the second the plasmon pole, which is the 
solution qp\{(jj) of the equation 

(«?)'- [1-(1-5'(9))V = 0. 

As in ID, the plasmon spectrum is gapless and makes an 
important contribution to low energy properties. 

In Eq. we need the imaginary part of Y{uj), which 
has a (S~function contribution at each of the poles of 
Y(uj)^ but also a regular part, and we find correspond- 
ingly Imy = y(P'^'') -I- r^Pi'^'^) + with the particle 
and plasmon contributions (for > 0) 
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y(plaB)(^) 



g[l-.gV^o^ 
2-52 
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and the regular part 



dq 



2v'^quj{l — g 



To proceed, we consider a bare interaction Vo{q) of the 
form 



27re" 



(1 



-2ijA 



which is the 2D Coulomb interaction in presence of a 
screening gate at distance A. Then Imy(a;) contains 
three energy scales: wq = 2i:vqv'^ = 4eF, = kw, where 
K = / £d is the two-dimensional inverse screening 

length, and = v/2A. For the bare 2D Coulomb in- 
teraction, i.e. A oo, only the two scales and Sl„ 
remain, that are related by Vs = 20^/^0 ■ For finite A it 
is convenient to introduce the ratio A = O^/^^a = 2kA. 

In 2D disordered systems Im Y diverges at low frequen- 
cies which leads to a divergence of J{t) for t ^ oo imply- 
ing a total suppression of !/(e) for e — > 0. The same is true 
for ID ballistic wires. In the 2D ballistic case considered 
here, Imy remains finite or even vanishes in presence of 
a gate as shown in Fig. ^ Accordingly, at T = 0, J{t) 
approaches a constant Jas for i — > oo, which is given by 



J.s = - — ImFM. 
Jo I" 



(13) 



It is important to note that the explicit result for lmY{uj) 
given above does not suffice to determine Jas explicitly, 
since the integral (|13|l also has contributions from fre- 
quencies that are not small compared to ep where the 
true parabolic dispersion must be used, which leads to 
a faster decay of lmY(uj) at high frequencies. However, 
this high energy behavior chiefly affects the quantity Jas , 




(olQ 



FIG. 1: The function linY(Lij) for an unscreened bare 
Coulomb potential (solid curve) and in the presence of a gate 
at distance A with kA — 5 (dashed curve) and kA = 1 (dot- 
ted curve). Y in units of tt/ojo = l/2uov'^, lo in units of 



while for large times the remaining part J'{t) in the de- 
composition 

J(t) = Jas + J'it) 

is determined by the low energy behavior of Y{uj). 
Accordingly, the spectral density splits into 

P{E) = SS{E) + P'{E) , 

with the regular part 

P'{E) := S9{E) [ ^ e'^* (e^'^*) - 1) , 
J 27r 



where S 



< 1 is a suppression factor. Hence, the 



(5(i?)-function form of the spectral density of a noninter- 
acting system partially survives in 2D ballistic electron 
systems. 

From Eq. (fTn|l the DOS now reads 

where the energy dependent part is given by 

iy'{e) ^uo /' ^ dE P'{E) . 
Jo 

Clearly, i'(e) has a nonzero value i^(0) = Si>q at e = 0, 
which depends on the integral property Jas and thus on 
the interaction strength as opposed to earlier predictions 
0,01 • Also iy'{e) contains the suppression factor S as 
a prefactor. Further, the slope and curvature of iy{e) at 
e = are not universal but depend on the interaction 
strength r^. 

To determine the value of Jas and S explicitly without 
artificial cutoffs, we either have to use the full quadratic 
dispersion or to approximate the less important regu- 
lar part of Y{uj). Since ^^''''^^(a;) approaches — l/4z/ow^ 
for small and large frequencies and does not deviate 
much from this limiting value in between, we replace 
it by that constant. In this approximation, which does 
not necessarily capture the true high frequency behav- 
ior and thus gives only an estimate of Jas when in- 
serted into Eq. H13() . the function lmY{uj) is of the form 
Imy(w) = ^y(Lu/D,f^, A). From this scaling form we find 
Jas = -rsC(A) with C(A) = 1/2 du y{u,X), which in- 



dicates that the factor S 



rapidly suppresses 



j'(O) with increasing rs- With decreasing gate distances 
the suppression factor approaches 1, which means that 
interaction effects become weaker. These features of the 
DOS are also apparent from Fig. 13 which displays results 
for various values of with and without gate. Without a 
gate the DOS is cusplike near e = 0, whereas in presence 
of a gate the slope at e = vanishes. This latter result 
is in accordance with the perturbative analysis • 

For finite temperatures and in the absence of a gate 
the denominator in Eq. ijHJl behaves singular at w = 0. 
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FIG. 2: Density of states for = 0.1, 0.3, 0.6, 1.0 (topmost to 
lowest graph). Inset: The correction to the DOS for rs — 0.6 
without a gate (solid curve) is reduced by a gate at distances 
A with kA — 5 (dashed curve) and kA — 1 (chain-dotted 
curve) . 



While \J{t)\ then exceeds |Jas|, it turns out that j/(0) 
always increases if the temperature is raised. This means 
that finite temperatures smear out the cusplike DOS and 
reduces the size of interaction effects near e = 0. 

Now, two remarks are in order. Strongly interacting 
semiconductor sys tems, where values for rg up to 40 have 
been observed [l5j, are usually disordered, but a ballistic 
regime is feasible |l6j. For large our analysis does 
not necessarily apply and, in particular, as was discussed 
above, we cannot obtain reliable values for Jag. However, 
if no scales other than ep and lire enter, Jas should be 
of the form Jas = ~'"sC('^s, A). Then S will vanish for 
Ts S> 1, provided limr,^oo ?'sC(^sj ''^) = 

Second, in the absence of a gate we have h'(e) « 
i'oS'[l + |e|/4eF + 0{e'^)]. This differs from previous per- 
turbative approaches giving z^(e) « z^o(l + |e|/4eF) jllLlia | 
independent of rg ■ The lack of a dependence on the in- 
teraction strength can be traced back to the fact that in 
these works Im Y{uj) is effectively replaced by a constant 
and cut off at the Fermi level. Then, no interaction scale 
remains and therefore the authors obtain a universal (in- 
teraction independent) result. From Fig. we see that 
this is not the case, but the dominant plasmon contribu- 
tion to lmY{uj) falls off on the scale £7^, which depends 
on the interaction strength. Indeed, if we turn off the in- 
teraction, the corrections to the nonintcracting DOS i^q 
vanish. 

Throughout this Letter, we have considered a pure 
2DEG without disorder. In weakly disordered 2DEGs 
with elastic mean time r electrons move diffusively on 
long time scales corresponding to excitation energies 
small compared to 1/t. Then, at T = the DOS drops 
down to zero for energies below 1/t where the behavior 
passes over to an Altshuler-Aronov ZB A [3-6] . Some au- 
thors have investigated the crossover from the diffusive 
to the quasiballistic limit perturbatively [lllll^ll7j |. The 



clean limit of these treatments does not yield our results 
for the reasons discussed above. 

Experimentally it is possible to realize very clean and 
strongly interacting 2DEGs (or 2D hole gases) [isl Il6| , 
but so far mainly transport measurements have been per- 
formed in the context of the metal-insulator transition. 
Tunneling experiments in samples with high have been 
suggested JlSj; but have not yet been performed, since 
they are more involved, because a counter electrode has 
to be added at a small distance without influencing the 
quality of the device. Nevertheless, tunneling experi- 
ments with a ballistic 2DEG seem to be feasible in the 
future. Interaction effects on the DOS in 2DEGs are also 
relevant for experiments with 2D-1D tunnel junctions 
[lij . because there not only the DOS of the ID wire is 
measured, but a convolution with the suppressed DOS of 
the 2DEG. 
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